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We identify the critical theory controlling the universal, low temperature, macroscopic properties of 
both quantum and classical ferromagnetic chains. The theory is the quantum mechanics of a single 
rotor. The mapping leads to an efficient method for computing scaling functions to high accuracy. 
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A number of recent papers have studied the finite 
temperature properties of ferromagnetic quantum spin 
chains. At low temperatures, macroscopic observables 
can be fully described by two dimensionful parameters 
which characterize the ground state. A convenient choice 
for these is the ground state magnetization density Mq, 
and the ground state spin stiffness ps. Then the macro- 
scopic properties of a quantum spin chain of length L, in 
the presence of an external magnetic field H, and at a 
temperature T are fully universal functions of the dimen- 
sionless ratios that can be formed out of these parame- 
ters. A convenient choice for these ratios is (recall that 
in d — 1 Ps, has dimensions of energy • length) 

r = psMo/T, h=H/T, q = ps/{LT). (1) 

Here we have used units in which h ^ — I and ab- 
sorbed a factor oigpB into H {pB the Bohr magneton). 
Thus, for instance, the temperature dependent magneti- 
zation density obeys jsj M — MQ^M{r,h,q) where 
is a universal function. This, and other, universal func- 
tions will depend upon the boundary conditions on the 
chain: we will focus, for simplicity, on periodic bound- 
ary conditions i.e. on spin rings. There is no method 
for computing in complete generality — in this paper 
we shall show how to compute $jv/ efficiently in a limit- 
ing case when the quantum ring can be described by an 
effective classical model. 

A convenient point for beginning discussion is the spin- 
wave expansion. This expansion is valid provided r ^ 1 
(provided h is not too small), and it is quite easy to use 
standard methods to determine the leading term. For a 
ring we find 

^Mir,h,q) = l-^J2 2 2 2 / 1 (2) 

r exp(47r^n^(j^/r + n) — 1 

where the sum is over all integers n. An interesting 
property of this expression emerges in the limit r — )■ oo, 
/i — > 0, but with 

g = rh = psMoH/T^, (3) 

fixed; then we find $m = 1 - (1/2^5) coth(^/2g). The 
implication of recent works 0-0] is that this limit is non- 



trivial at each order in the spin- wave expansion, and that 
the resulting series has in fact properties of the classi- 
cal ferromagnetic ring. Thus we may define the classical 
scaling function 0^/ by 

0Af(5,9) = lim «'a/(?',5/^,9)- (4) 

r — *oo 

Classical behavior emerges in this limit because the fer- 
romagnetic correlation length becomes larger than the de 
Brogiie wavelength of the spin waves. 

One possible approach to the computation of the scal- 
ing function (pm is to compute the magnetization of a 
nearest-neighbor, classical ferromagnetic chain, whose 
statistical mechanical properties were computed some 
time ago [^-^. The scaling limit of classical solution 
was studied in recent work |^,^, and led e.g. to the re- 
sult (t>Mig, 0) = |5 — -^g^ + 0{g^) — this result means 
that the usual linear susceptibility dM/dH diverges as 
T^^ and that the third order non-linear susceptibility 
d'^M/dH^ diverges as T^^. However, the computations 
required to achieve this limited result were quite compli- 
cated. In this paper we shall develop an efficient method 
to computing the complete function (f)M{g,q) to essen- 
tially arbitrary accuracy. This will be done by a pre- 
cise identification of the critical theory controlling this 
crossover function. 

We begin by considering the partition function of a 
classical Heisenberg model in a uniform magnetic field: 

where the exchange constants J{i) > 0, J(i) — J{N — i), 
the Ui are unit 3-component vectors which are integrated 
over. The universal critical theory emerges when we take 
the continuum limit of this action. With a lattice spac- 
ing a, the continuum limit will be characterized by the 
values Mq = 1/a, L = Na, and ps = aYl^i^ J{i). Our 
results apply to all systems in which the summation in 
the definition of ps is convergent — for J(i) ^ i^'^ this is 
the case for q > 3. The continuum field theory which 
emerges by this method gives the partition function Zc: 
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where the integral is now a functional integral over unit 
vector fields n(a;) satisfying n(0) = n(L). A key property 
of Zc is that it is a finite field theory, free of ultraviolet di- 
vergences. This becomes clear when wc re-intcrprct Zc as 
the imaginary "time" {x) Feynman path integral for the 
quantum mechanics of a single particle with co-ordinate 
11(3;) restricted to lie on a unit sphere: no discretization 
of time is required to define this problem. As a result, 
all observables are universal functions of the couplings in 
Zc, and it is useful to now transform to dimensionless 
variables. We rescale spatial co-ordinates y = Tx/pa, 
and obtain Zc. 



Subsequent computations are best carried out using 
the Hamiltonian, H, of the quantum particle described 
by ^c: 
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(5) 



This describes a single quantum rotor with unit moment 
of inertia, angular momentum operator L (which obeys 
the usual commutation relations [L^^L^] = ica^-yL^), in 
the presence of a "gravitational" field g. There is no 
need to consider the radial motion as the length of n is 
constrained to unity. The logarithm of Zc equals the free 
energy of the quantum system W at a "temperature" q; in 
the original spin ring, q is the ratio of correlation length 
at = to length of the system. For other boundary 
conditions, Zc will be given by appropriate propagators 

ofn. 

We now consider eigenvalue equation Hip = Eip. As 
H commutes with L^, eigenstates are divided to sub- 
spaces of azimuthal quantum number m = ... — 3, 
—2,-1,0,1,2,3,.... In spherical co-ordinates {9,^), ip 
is given by e""*'w(^). At g = u{9) is given by the asso- 
ciated Legendre Polynomials F,'"''(cos0) and eigenstates 
of H are spherical harmonic states \l,m) with I > \m\. 
The matrix elements of 7i in this basis are 



I'+i 



4P 



1 



+1 



;'2 



41'' 



1 



(6) 



Notice that H is tri-diagonal in each m subspace: this 
makes numerical diagonalization of H. quite straight- 
forward. The eigenvalues of Hamiltonian are given by 
Em,n{g), and n = 0,1,2,... is the number of nodes of 
function u{9). We also generated a power series expan- 
sion in g for the ground state energy Eofi{g) using a sym- 
bolic manipulation program (Mathematica) ; this leads to 



the magnetization scaling function for the infinite ferro- 
magnetic ring (/)m(<?,0) = —dEofi{g)/dg, for which we 
find 
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All coefficients are rational numbers. 

In the complementary large g limit, the particle spends 
most of its time near the 'north pole', and in its vicinity 
it experiences a harmonic oscillator potential well. It 
therefore pays to work now in the basis states of this 
harmonic oscillator and thereby generate a perturbation 
expansion valid for large g. Parametrizing E = —g + 
y/ge and u{e) = (6i/sin(6'))V2/(gi/46i) ; this gives us an 
eigenvalue equation for /: (ho + h\)f{z) = ef{z) with 
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where z = g^^'^9. Notice that fto describes a two- 
dimensional harmonic oscillator in radial co-ordinates. 

Its eigenstates n > have energy Eq = 2n + |to| + 1 
and are represented by generalized Laguerre polynomi- 
als 0l™lL[^™'|^|(z2)e-^'/2_ Further, notice that hi can 
be expanded as a series in positive integer powers of 
z^, with all terms being small for large g. The matrix 
elements of hi in the |n) basis can be determined by 



repeated use of the identity z \n) = (2n 



l)|n) 



-^Jn{n + \m\)\n - 1) - ^{n + l)(n + \m\ + l)|rH- 1). It 
now remains to diagonalize h in the \n) basis, which can 
be done order by order in g~^/'^ by Mathematica. Such 
a procedure was used to generate an expansion for the 
ground state energy, i?o,o(.9) and hence for (puig, 0): 



<f>M{g,0) 
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(8) 



Unlike the small g expansion, which has a finite radius of 
convergence, the large g expansion is only asymptotic. In 
particular, the large g limit loses topological information 
associated with tunneling paths which traverse the south 
pole — such paths will lead to 'instanton' contributions 
which are exponentially small for large g. In Table 1 we 
give the higher order coefficients of these expansions. 
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For the quantum ferromagnetic Heisenberg ring with 
spin 1/2: 
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we can calculate the magnetization at A'^ = oo limit for 
given temperature and magnetic field using thermody- 
namic Bethe ansatz equations 1^. The magnetization 
still obeys the same limiting scaling function. The stiff- 
ness constant ps is Ja/4 and Mq is l/2a. In Fig. 2 we 
compare the magnetization as a function oi g — JH/8T^. 
As T goes down, the magnetization divided by Mq ap- 
proaches the scaling function (/)Af((7,0), supporting the 
conclusion that the quantum ferromagnetic Heisenberg 
ring has the same magnetic scaling function with classi- 
cal one. 

In our scaling theory, the linear susceptibility dM / dH 
of quantum ferromagnet diverges as at low temper- 
ature. Schlottmann ||^ proposed the divergence of the 
type T~^/ ln( J/T) using the numerical analysis of ther- 
modynamic Bethe ansatz equations. But later this was 
negated by the detailed investigations of thermodynamic 
Bethe ansatz equations pOj . 

Next we consider the function ^jv/Cff, q) at q > 0, which 
is important for the analysis of short rings. This is rep- 
resented by the eigenvalues Em.n{g) of Hamiltonian (||): 
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Em En Ka,n{9) exp(-£;,„,„ (ff) /g) 

Em E„ exp{-£;„,„(c/)/(7) 



(10) 



For (7 ^ 1, is dominated by the ground state m = 
n = 0. The energy gap to the second lowest eigenvalue 
at TO = ±1, n = is more than 1. Then, deviations 
from (Pm{9,0) are exponentially small: 



^Mig,q)^^M{9,0) + Oie-'/''), <7«1 



(11) 



In Table 2 we give the result of numerical calculation of 
4>M{g, q)- We calculate Em,n{g) and E'^ ^{g) numerically 
by diagonalizing the tridiagonal matrices (^). Terms at 
very big n or to are not necessary because their contri- 
butions are exponentially small. 

Finally, we note that the behavior of the scaling func- 
tions is also simple in the limit g ^ 1. The problem is 
now equivalent to a single classical rotor: 

(l>M{g, q) = coth , - = — - — . (12) 

q 9 q T 

This means that the system behaves as one big spin MqL 
if the correlation length ps/T is much longer than the 
system size L. 
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TABLE L 
g expansion ^ a„(/ 



Expansion coefficients of 
^""^ and large 



m(5,0) 
g expansion 1 



for small 



n 




bn 
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0.66666666666666667 


-0.0078125 
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-0.005859375 


3 


0.265255731922398589 


-0.004852294921875 
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-0.243362205238748449 


-0.0045318603515625 
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0.235324701717019961 


-0.0047218799591064453 
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-0.234382743316652222 


-0.0054293274879455566 
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-0.0068305558525025845 
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-0.24474146816049874 


-0.0093398639000952244 
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-0.0138054155504505616 
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-0.02195840459989995 
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-0.037433421774196063 
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-0.296386895761683996 


-0.068149719593456837 
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-0.132063536806409254 
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-0.33555167675962067 


-0.271575442615252266 
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0.35866981723790861 
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-0.38433633684212508 




17 


0.41273495146936857 
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-0.44407985807794498 
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0.47861575559844195 
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TABLE 11. Values of scaling function 4>M(g, q) 
for q = 0.0,0.5,1.0,1.5,2.0 
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FIG. 1. The scaling function <j!>Af((?, 0). Line a is expansion 
(^) up to g^'^ and Line b is expansion up to g^'^. Line c is the 
result of asymptotic expansion (0) up to g^^^^^ from g = oo. 




FIG. 2. Magetization versus g = JH/8T^ for spin-half and 
infinite-length ferromagnetic Heisenberg chain with nearest 
neighbor exchange (^. As temperature goes down, the line 
approaches to the theoretical line 4>M(g, 0). 




FIG. 3. Scaling function <j)M{g,q) for various values of q. 
Solid line is for 5 = 0, dashed line is for q — 0.5, dotted line 
is for g = 1.0 and dashed chain line is for q — 1.5. 
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